Angular momentum changing transitions in 
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ABSTRACT 

Collisions between electrically charged particles and neutral atoms are cen- 
tral for understanding the dynamics of neutral gases and plasmas in a variety of 
physical situations of terrestrial and astronomical interest. Specifically, redistri- 
bution of angular momentum states within the degenerate shell of highly excited 
Rydberg atoms occurs efficiently in distant collisions with ions. This process is 
crucial in establishing the validity of the local thermal equilibrium assumption 
and may also play a role in determining a precise ionization fraction in primor- 
j dial recombination. We provide an accurate expression for the non-perturbative 

rate coefficient of collisions between protons and H(n£) ending in a final state 
}l{ni'), with n being the principal quantum number and i, £' the initial and final 
angular momentum quantum numbers respectively. The validity of this result 
is confirmed by results of classical trajectory Monte Carlo simulations. Previ- 
ous results, obtained by Pengelly and Seaton only for dipole-allowed transitions 
£ — )■ £± 1, overestimate the ^-changing coUisional rate coefficients approximately 
by a factor of six, and the physical origin of this overestimation is discussed. 

Subject headings: atomic processes — stellar astrophysics — photosphere — early 
universe 



1. Introduction 

Collisions between electrons or protons and Rydberg atoms modify the atomic level pop- 
ulations, enhancing or suppressing statistical equilibrium, and upsetting selection rules and 
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decay channels established by purely radiative transitions. Electron-H(n) collisions are more 
efficient in transferring energy than p-H(n) collisions, while the latter are far more effective 
in mixing £ states. CoUisional capture of electrons into highly-excited states is the primary 
process by which atoms form in cold and ultracold plasmas, and by which anti -hydrogen 
atoms are created in non-neutral magnetized antiproton and positrons plasmas (IGabrielse 
20051 ). During this stage, high-£ states are preferentially populated, and subsequent collisions 
with electrons (and ions) populate low-£ states that have much faster radiative de-excitation 
rates. 



Early theoretical results for proton -hydrogen atom collis i ons we re obtained by lDemkov et al. 



( I1974J ) for low quantum numbers, by iPercival fc Richards! ( 119771 ) using a classical model, 
and by lAbrined (Il966[ ) employing classical trajectory Monte-Carlo (CTMC) simulations. 
Measurements of low ^-m ixing rates were performed for beam gas Na^-Na(?T,£) collisions 
( Sun fc MacAdamd Il993l ) and theoretical comparison was provided via Floquet analysis 
( ICavagnerd I1995I ) . More recently, electron ^-changing collisions were shown to be responsi- 
ble for transferring magneto-optically trapped rubidium Rydberg atoms into high angular 
momentum states (IDutta et al.ll200ll ). In this case, electron ^-mixing efficiency is boosted 
by the presence of the trapping potential allowing for frequent collisions with the Rydberg 
atoms. The £- changing processes are also important in Zero Electron Kinetic Energy (ZEKE) 
spectroscopy (ISchlag &: Levind 119971 ). 



In a more strict astrophysical setting, coUisional physics, both n-changing and f-changing, 
is critical to understand the validity of the assumptions on which the determination of stellar 
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equilibrium applicable in the stellar interior - where the mean free paths for photons, elec- 
trons, and ions is rather small with respect to their average distance - no longer holds near 
the surface of the star. In an intermediate regime, the photons escape with a large mean 
free path, while electrons and ions still maintain small mean free paths for collisions. This 
is the origin of a peculiar status in which, although the photon distribution departs from 
the equilibrium black-body distribution, the gas particles still achieve energy distribution 
characteristic of thermodynamic equilibrium. This allows for the evaluation of gas particle 
population via the knowledge of the local temperature alone, known as Local Thermody- 
namic Equilibrium (LTE). 

However, in the upper part of the stellar atmosphere, the chromosphere, the particle 
densities drop so much that even collisions among the gas particles and photons are not 
enough for maintaining LTE. This leads to a breakdown of LTE, with the atomic popula- 
tions and the ion and electron densities no longer determined by the Boltzmann and Saha 
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formulas, and the detailed knowledge of atomi c level populati on for each species, including 
angular momentum states, is instead required ( ISampsonlll977 ). The validity and the break- 



down of LTE depend upon the strength of the atomic transitions. For transitions having 
optical depths larger than the optical depth in the surrounding continuum, the thickness of 
the corresponding region is relatively small, allowing one to meaningfully use the average 
values for the entire region. This is not true if the line-forming region is located in the 
chromosphere, where the condition of LTE has to be relaxed. Collisio nal excitations are als o 
important for the determination of the primordial helium abundance (ILuridiana et al.ll2003l) . 



the understanding of line formation for early-type stars (IPrzvbilla fc But 
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Moreover, the redistribution of ^-state populations in B.{n£) in high redshift universe 
can change the spectrum of the primordial recombination at low frequencies, leading to emis- 
sion of successively lower-frequency photons and spectral distortion. The new surveys of the 
cosmic background background, such as the 7-year data integration of the WMAP satellite, 
and the forthcoming data from the Planck Surveyor, promise higher-precision determina- 
tion of the cosmological parameters, in particular the spectral index of scalar perturbations 
and the baryon content of the Universe. These two parameters are directly affected by 
modifications to the cosmological recombination models. Mixing among high-£ states, for 
lower redshifts, leads to an increase in recombination to the grou nd state, while n-changiri g 
collisions suppress the emission of recombination epoch photons (IChluba et al.l 120071 . l2010l ). 



In this work, we provide analytical expressions for ^-mixing rate coefficient in collisions 
between a Rydberg hydrogen atom and a proton, represented by the process: 

H(n, £) + H+ ^ H(n, £') + H+. 



The theoretical model can treat the more general case of collisions with an ionic projectile 
of mass M and charge Z. A semiclassical perturbation theory for £ — > £ ± 1 collisions 
( jPengelly fc Seatonlll964l ) is divergent for both small and large impact parameters, requiring 
the introduction of an ad-hoc radial cut-off to regularize the results. The model proposed in 
this paper is non-perturbative and the only assumption made is that the projectile moves on a 
straight line. This is well justified by the fact that collisions with large impact par ameter have 
large probability for ^-changing, as first shown in lVrinceanu fc Flanneryl ( l2001al ). Moreover, 
our model predicts that the rates decrease roughly as |A£|~^. 
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2. Angular momentum changing transitions 

In our model for collisions between charged particles and Rydberg atoms, we make three 
basic assumptions. First, we assume that collisions occur on timescales which are much 
longer than the orbital motion of the Rydberg electron. This implies that the Rydberg 
electron travels many orbitals during one collision in order to appreciably change its angular 
momentum. Second, we assume that the colliding particle creates a weak electric field, 
yet large enough with respect to the Stark electric field required to lift the hydrogenic 
degeneracy in each Rydberg manifold of principal quantum number n. Third, for pure 
angular momentum mixing without energy exchange, it is necessary that collisions occur at 
large impact parameters. Within the above mentioned approximation scheme, the physical 
picture for ^-changing collisions is the following. For collisions at large impact parameters, 
it is reasonable to assume that the charged projectile move s mostly undisturbed, along a 



classical straight line trajectory, creating a weak electric field (IVrinceanu fc Flanneryll2001al ). 
Under the influence of this electric field, the Rydberg electron slowly precesses mixing its 
angular momentum state inside the degenerate shell. This collision is so weak that the 
angular momentum changes without any exchange of energy. Because of the long range 
nature of the Coulomb interaction, the cross section for this process can be very large, and 
indeed we will find that it is logarithmically diverging for AC. = ±1 transitions, requiring a 
semiclassical treatment to cure this divergence. Based on geometric arguments, it is expected 
that the cross section scales as a ~ vrn^aQ, where Oq is the Bohr radius. The results introduced 
in this paper are valid for n > 10 (semiclassical approximation), for projectile velocities 
smaller than the orbital velocity of the Rydberg electron, which implies an upper bound on 
ny/T < 2.4 X lO'^K^^'^, and for densities greater than Ncrit as defined by Eq. f[T2|) . 



2.1. Straight-line trajectory time dependent probability 

Within the above mentioned approximations, the time dependent Schrodinger's equation 
can be analytically solved for the evolution of the states withi n a Rydberg shell, subject to th e 



electric potential created by the passing charged projectile ( IVrinceanu &: Flanneryll2001al ). 
The key step in finding the exact solution is the observation that the angular momentum 
and Runge-Lenz operators generate an SO (4) symmetry group, which in turn decomposes 
into the direct product of two rotation groups SO(3)(g)SO(3). 

The cross section for ^-changing collisions within an energy shell with principal quantum 
number n can be written as an integral over the probability of the impact parameter as 

POD 

atl, = 2n Pi^bdb, (1) 
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where the probabihty for making £ — > transitions has the form, 



n-l 



L=\l'-t\ 



(2) 

where {■ ■ ■ } is a six-j symbol, L is the vector sum of ^ and is the ultraspherical 

polynomial, and j = (n — l)/2. Th e rotation angle x is defined as (jVrinceanu fc Flannery 
2001al : iFlannery fc Vrinceanull2003l ) 



cosx 



1 + cos(A<I>a/1 + a^) 



(3) 



l + a2 

where A$ is the azimuthal angle swept by the projectile, such that in a complete collision 
A$ = TT, and a is the scattering parameter 



a 



3 nh 

— 7 

2 mgf 



(4) 



where mg is the electron mass and v the initial velocity of the colliding projectile. The scat- 
tering parameter is directly related to the ratio between the maximum angular momentum 
allowed by the Rydberg electron and the initial angular momentum of the charged projec- 
tile. A c omplete derivation of the abo ve expressions can be found in IVrinceanu fc Flannery 
f |2nnia[ ): [Plannerv fc VrinceamJ (120031 ). 



The probability (|2]) depends on the impact parameter b and the projectile velocity v 
through a, such that the cross section in Equation ([1]) can be rewritten as 

2 



cr 
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where the velocity independent integral factor is determined by the initial and final 

states as: 



An) 



pin) , X da 



(5) 



Since the collision probability scales as P^^i oc a^'^"^'', and the rotation angle x ~ 2a 
for small a, the integral factor I^^i^i-, and therefore the cross section, diverges for collisions 
in which |A£| = 1, i. e. the dipole allowed transitions. The logarithmic singulari ty in 
the cross section for large impact parameters is well known (jPengelly fc Seatorull964j ) and 
is a refiectior i of the impact paramete r and Born approximation. This difficulty has been 
addressed in iPengelly fc SeatonI ( 119641 ) by considering either many-body collective effects, 
such as the plasma screening effects or Debye shielding, or the Rydberg natural linewidth, 
that limits the duration of ^-mixing collisions. 
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2.2. Semiclassical probability 

The non-perturbative probability P^^^i in ([2]) can be computed with the desired numer- 
ical precision for a range of quantum numbers, typically n < 60, before becoming unstable 
due to the oscillatory behavior of the ultraspherical polynomials. In order to carry out calcu- 
lations efficiently and analytically beyond intermediate Rydberg principal quantum numbers, 
we use the asymptotic limiting forms of the six-j coefficient and the ultraspherical polyno- 
mials, CiC\ and obtain analytical expressions for the rate coefficients in the limit of large 
quantum numbers. A derivation of these results is detailed in the Appendix. The probability 
for ^-mixing transitions in this semiclassical (SC) limit can be written as 



P''^{i/n,i'/n,x) 



2f 



Tchn'^ sin x 



, if I sinxl < I sin(?7i — ri2)\ 
KiB/A) I . I I . , , ,1 



A- , if I smxl < I sm(r7i +772)1 

V-D 



where K is the complete elliptic integral, A = sin^x — sin^(?7i — 772)? B = sin^(?7i + 
p-}) — sin^(r7i — r]9,), cos 771 = £/n, and cos 772 = (^'/n. The same result was obtained by 



Vrinceanu fc Flanneryl ( 120001 ). based on calculating the overlap between volumes in the clas- 



sical phase space. 

Notice that the scattering probability is zero for sinx < | sin (771 — 772)!, which implies 
a sharp cutoff for and consequently for the minimum allowed value of a for which the 
transition is possible, amin- The integral 



POO 



da 
a^ 



can be obtained in the semiclassical approximation and is divergence-free even for \Ai\ = 1 
transitions. Moreover, the calculation of the integral simplifies for small angular momentum 
transfers Ai. In this case. Equation ([7]) is obtained by writing cos 7/1 = i/n and cos 7/2 = 
{i + Ai)/n and calculating the Laurent series expansion in powers of Ai. In this process, 
we take advantage of the fact that the elliptic function at origin is K{0) = tt/2 and that 
sinx ^ 2a for small a. This simplified expression for transition probability can be easily 
integrated to obtain a non-perturbative analytical expression for the corresponding cross 
sections and rate coefficients, for any Ai. 

In Figure [H the collision probability P^^^, oc a^'^^^'l is plotted versus the impact param- 
eter b, scaled in units of the Rydberg size a„ = ra^ao, for different approximation schemes and 
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a specific dipole (A£ = 1) tr ansition (n, £ = 40 , 36 — )■ n,i' = 40,35). We can compare the 



Born approximation used in iPengelly fc Seatonl (119641 ). the expression in ([2]), and its semi- 



classical approximation given by Eq. El The result from a direct calculation of Equation ([2]) 
contains oscillations around the base line formed by semiclassical approximation Equation 
([6]), for most of the impact parameter range. Figure [1] shows that the probability for an- 
gular momentum mixing grows linearly with the impact parameter within the semiclassical 
approximation. This is a counterintuitive conclusion, since one would expect that the ability 
of a projectile to induce changes in a target atom would decrease as the impact parameter 
increases. Due to long range forces and zero energy exchange, angular momentum chang- 
ing collision are most effective at large impact parameters, even for separation of several 
hundreds Rydberg atom radii. The semiclassical approximation falls abruptly to zero at a 
critical impact parameter given by the condition | sinxl < sin(?7i — 772) in Equation (E]). The 
quantum expression in Equation ([2]) has a gradual decrease to zero, eventually merging with 
the Born result for sufficiently large b, as or a^. This is the source of the logarithmic 
singularity discussed for Equation ([5]). 



The transition probability in iPengelly fc Seatonl (jl964j ) is inversely proportional to 6 



1.2 



posing two difficulties: at small b the probability diverges and therefore it cannot be prop- 
erly normalized, and at large b it leads to a divergent cross section. Therefore two cut-off 
parameters are required. A short cut-off Rc is introduced at the point where the transition 
probability is 0.5, under the assumption that for b < the probability has fast oscillations 
that average to 0.5, replacing it by a constant. The more serious divergence at large impact 
parameters is cured by the second cut-off, that takes into account either the finite density 
of the medium (plasma) where such collisions occur or the finite time of these collision due 
to the radiative lifetime of the Rydberg atom. Figure [1] shows that the Born cross sections 
overestimate the quantal and semiclassical cross sections. Our quantum calculation leads to 
divergent cross sections only for Ai = ±1, while the semiclassical approximation provides 
finite results for any transitions, without the need for a cut-off parameter. 

Further simplifications, accurate for most cases of interest, occur if the scaled angular 
momentum i/n in the semiclassical formula ([6]) is assumed to be a continuous parameter with 
values ranging between and 1, and the the probability is expanded in powers of \A£\/n. In 
this case the integral factor can be calculated term by term, starting with l/\Ai/n\^, to get 



1£< /2[l-(£</n)V3 , [1 - 3(£</n) V(3^</^) , ^ 



n e { \i'/n-e/n\^ \e'/n-£/n\'^ \\e'/n-£/n\ 

where = min (£,£'). 



(7) 
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In Figure[2], we compare the integral factor for transitions inside the n = 40 shell, starting 
from I = 8, 20, and 32 and going to all possible final angular momenta. Three methods 
are used in this calculation: integration of the quantum formula direct integration of 
the semiclassical result ([H]), and the asymptotic expansion ([7]), dominating small angular 
momentum transfers. The simple expression ([7]) is quite effective even for larger transfers. 
The next term in the expansion, proportional to 1/Ai, is easily obtained, but it does not 
significantly improve the agreement with quantum calculations. 



2.3. Rate CoefRcients 

The rate coefficient for ^-changing transitions in contact with ions having a Maxwell- 
Boltzmann distribution /mb, is 



V a}^^, tMB{v)dv = 1 ^ 1 y 



2kBT 

where M is the reduced mass of the ion-atom system, and T the temperature of the ions, 
hereafter assumed to be expressed in Kelvin. 

By using the simplified semiclassical expression ([7]), the coUisional rate coefficient for 
ion induced transitions is 



TT , fW n^n\£ + f)-£l{i + £' + 2\Ai\)] 

Qne^ni'K-^ ) — -JA / — — 5—5- \ T?l — Tf; ^ 



2 mle'^ \l h?kBT \ (£+l/2)|A£|3 
The rate coefficient in cm^/s is then 

W n^n\i + £') - ilii + i' + 2\Ae\)] 




,..,^„„^ 1.294 X 10- (/+l/2)|Ay 

As expected, this rate scales as with the principal quantum number and rapidly drops 
with the increase in angular m omentum c hange , as 1/|A£|'^, similarly to the A£ > 1 fitting 



formula in Equation (4.112) of iBeigmanI (119951 ). The maximum rates are obtained for the 



dipole allowed transitions, £ — )■ £ ± 1. 



It is instructive to compare our two-body rate coefficient with the results of iPengelly &: Seaton 



(119641 ). which only accounts for dipole-allowed transitions. 



nPS(T\ mee^ 2 fM' PS 

Inl [T) = W - -— W-r^T-^ ^ \ — ^nl X 



, / N / 1 — 7 , k^rrie , TrUe , T 
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where 7 = 0.577216 is the Euler constant, the coefficient is D^f 
Nf. is the electron density expressed in cm~'^. 



6n2(n2-£2_£_i)^ and 



Our rate coefficient in Equation ([H]), summed over the two dipole transitions is 



(lnl{T) 




where the coefficient in this case is Dnt = Qn^iin? — i'^ — 1/4^). The ratio qn£{T) / q^f (T) 
becomes a slowly varying function of n and i. For example, for (T = 3, 000 K, = 300 
cm~^) this ratio is 6.98 for {n = 50, i = 48), and 6.20 for (n = 100, £ = 98). The reason 
for which the Pengelly and Seaton formula overestimates the angular momentum changing 
rates lies in the Born approximation. 

In order to compar e the ^-changing c o llision s with radiative processes, the quantity 
QniTnt was iutroduccd in lPengelly fc SeatonI (119641 ). where g„£ = qni^ni+i + Qni^ni-i and r„^ 



is the radiative lifetime, which is approximately given by r„^ ~ lO'^^n^i"^ s. This quantity 
is the inverse of a critical density above which ^-changing collisions are faster than radiative 
decay from level ni, which in our case is expressed through the relationship (in cm~^) 



crit 



nl) 



73 




(n/40)- 



M (^Vn2) (1 - £Vn2 



:i2) 



If scaled d own by a factor of appro ximately six. Equation ( !T2|) agrees with previous 
Born results in lPengelly fc SeatonI ( 1l964j ). 



We may also compare the rates for angular mo mentum mixing with other processes 
in plasmas of astrophysical relevance, as discussed in iDalgarnd (jl983[ ). Firstly, the angular 
momentum changing rate in collisions with electrons is smaller than the rate in collision with 
protons by factor of (mp/me)^/^ ~ 43. The energy changing rate due to proton collision is 
smaller than the angular momentum changing rate by [An/n)^. This factor comes from the 
assumption of straight line trajectory, which implies a transition probability proportional to 
the square of the dipole matrix element. The matrix elements for in-shell transitions are 
greater than those for inter-shell transition by approximately the same factor. 

The electron excitation rate coefficient depends as ~ j^5/6^i4/3 t he temperature T 
and n, while the de-excitation rate coefficient depends as ~ T~^/^n^/^ ( Pohl et al.l I2OO8I ). 
Therefore for low enough temperatures, these processes will be dominated by the proton 
induced ^-mixing. The radiative relaxation rate, represented by the inverse of the radiative 



life time r„£ ^ 10" 



10^3^2 



s, does not depend on temperature or density. A critical electron (or 



proton) density can be defined such that the proton induced ^-mixing rate is faster than the 
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radiative rate ( iPengelly fc Seatonlll964j ). A near-resonant charge transfer is most efficient at 
projectile velocities matchin g the velocity of Rydb erg electrons and at impact parameters 
in the range of (5 — 10)?t,^ (jSmith fc Chupkal Il995l ). This process is therefore unlikely for 
^-changing collisions which occur with large probability at non-velocity matching conditions 
and much larger impact parameters, see Fig. 1. Our CTMC simulations con&m that, in the 
range of parameters considered here, the number of trajectories resulting in charge transfer 
are a small fraction (less than 10%) of the total number of ^-changing trajectories. We 
compare rates for these processes at different temperatures in Table [H showing that the 
proton induced ^-changing is the fastest process over a wide range of parameters. 



3. Monte Carlo Simulations 



In order to conffim the range of validity of our semiclassical results in Sec. 12. 3[ we 
perform numerical simulations of proton (projectile) and H(n£) (target) collisional ^-mixing 
and calculate the rate coefficients for these collisions with the CTMC techniques. The 
earliest applic ation of CTMC technique to proton - hydrogen atom collision is reported in 



Abrined (Il966[ ). who calculated the classical charge transfer and ionization cross sections inp- 
H(n = 1) collision, with the initial ground state kinetic energy sampled from a microcanonical 
distribution. We similarly use a microcanonical population for the Rydberg atoms with 
generic principal quantum number n and angular momentum quantum number i. The 
velocities of the projectile ion and the center-of-mass of the Rydberg atoms are sampled from 
a Boltzmann-Maxwell distribution at temperature T. The collision is started at t = —tmax 
and ends at t = tmax where tmax = V^/'^i2, where r] = 4 is chosen in our simulations, 
and Vi2 is the relative velocity between the projectile and Rydberg atom. We perform the 
simulations at temperatures large enough such that the propagation time is manageable, 
but small enough such that energetic collisions could be neglected. All 18 conjugate degrees 
of freedom are propagated in time with a time-adaptive implicit Runge-Kutta integration 



Table 1: Rates in s ^ for processes occurring in a plasma with number density 10^ cm ^ and 
various temperatures 



Process 


10,000 K 


1,000 K 


25 K 


H"*" angular momentum changing 


6362 


20119 


127241 


Electron de-excitation 


5725 


840 


38.9 


Electron excitation 


14.7 


21.6 


40 


Radiative relaxation 


2442 


2442 


2442 
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method of forth order. The selection of trajectories for the calculation of rate coefficients is 
achieved by ensuring conservation of the total energy and by maintaining the on-shell energy, 
i. e. \n' — n\ < 0.5 where n' is the effective quantum number n' = l/\/—2E' for both the 
target, or the charge transferred atom, at the end of the collision. If the angular momentum 
of the direct or charge transferred Rydberg atom at the end of the collision is in the interval 
{£',£' + 1), then the trajectory is counted as a collision with final angular momentum i'. 

Care must be exercised in the choice of impact parameter distribution, considering 
the large range of b values which are required for effective angular momentum mixing, as 
illustrated in Figure [TJ According to the model, the probability for a change of Ai grows 
linearly with 6 up to a maximum impact parameter bmax = 3n^e/ (v| A£|), and is zero for b > 
bmax- There are two options to address this problem; one is using the importance sampling, 
i. e. to try to get a distribution that resembles the integrand, using more points where the 
integrand is large, or the option of using stratified sampling, in which the integration domain 
is partitioned into non-equal slices, each having uniform probability distribution. For the 
stratified sampling at a given Vi2, the b space is divided into intervals Di, D2, ■ ■ ■ , D^, where 
each segment is defined by 3n^e/{k + 1/2) < bvi2 < 3n^e/{k — 1/2). If our model is correct, 
a collision with impact parameter in a segment Dk can lead only to collisions for which 
\Ai\ < k. 

We have performed simulations for 3 cases: (a) n = 20, i = 4, T = 800, 000 K, (b) 
n = 20, i = 4, T = 400, 000 K, and (c) n = AO, i = 8, T = 150, 000 K. For cases (a) and 
(b), we run simulations for segments Di to Diq, and Di to Dig for case (c), with 40,000 
trajectories within each segment. Figure |3] shows the results of the simulations. 

In order to represent the results from these cases on the same graph, we scale the rates 
by 6n-^7rM/ (2/cbT), and the angular momenta scaled hy n, X = i/n, X' = i'/n. By scaling 
with the same factor above the rate coefficient. Equation can be written, in atomic units, 
as 

K(A.V) = |/S„ = ^ 

which does not depend on temperature, and depends on n, i and i' only through A and A'. 

Figure [3] shows that the ^-changing rate is divergent at i'/n = 0.2, because this corre- 
sponds to \A£\ — )■ 0, and decreases as 1/|A£|^ for larger Ai. The agreement between the 
simulated cases and the prediction of the simple formula (Q confirms the validity of the 
semiclassical model. The two curves are slightly different at small Ai because they actually 
represent the scaled rate (I13p integrated over the finite bins used in the simulations. The 
value of the scaled bin Si' /n is different for case (c). We have also verified through CTMC 
simulations that ^-changing rates for electron- Rydberg B.{ni) are smaller with respect to the 



A 



lAAl 



+ 



3X1 



2A^|AA|2 



(13) 
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proton ones by a ratio ^Jmyjme — 40, as expected from Equation ([8]) 



4. Conclusions 

Analytical expressions for the rate coefficient of ^-changing collisions between protons 
and H(n£) Rydberg atoms have been derived using a non-perturbative approach free of 
divergences, plaguing previous results for the dipole-allowed ^-mixing collisions. The results 
have been compared to CTMC simulations over a range of temperatur es of astrophysical 



interest. The dipole-allowed ^-changing collision coefficient evaluated in lPengelly &: Seaton 



(119641 ) overestimates the corresponding rates by about an order of magnitude. 
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Appendix: Derivation of Equation (6) 

In this section, the semiclassical expression (jH]) is derived from Equation ([2D in the limit 
of n — )■ oo, and the same scaled angular momenta ijn and i' jn. 

By using Euler-Maclaurin formula, the summation of Equation (|2]) is approximated by 
an integration such that 



2 

hm = 2i'n f{[ [ \ \ Hlix) d{^)+o(l] (14) 



I J J J J J \n 



where Hjl is th e generalized character associ ated with the irreducible representation of the 



rotation group ( IVrinceanu &: Flanneryll2001bl ). defined by: 



\2j + l){2j-L)\ ^^^^_^^LMi+L) 



HM = L\J (2J + L + 1)! ' ^^^^ 



in terms of ultraspherical polynomials Cn""*. In this section j = {n + l)/2. 



The classical limit of the 6-j symbol was first discussed by IWignerl ( 119591 ) . based on the 



physical interpretation that the square of the 6-j symbol gives the quantum probability of 
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coupling 3 angular momenta ji,j2 and to the total angular momentum j. The correspond- 
ing classical probability is then linked to the volume of a tetrahedron that has ji, j2 and 
js as edges joining at a vertex, and j, ji2 and J23 as the other edges, opposing that vertex. 
Here ju and J23 are the intermediate coupling angular momenta. This provides only an esti- 
mate, and the agreemen t is quantitative on l y on a verage over neighboring angular momenta. 
Heuristic arguments led iPonzano fc Reggd (119681 ) to an improved formula which takes into 
account the oscillations inside the classically allowed region, but which fails in the vicinity 
of the turning points. A rigorous derivation of these results and an uniform approximation 
valid over a larg e range of values of the ang ular momentum, including classically forbidden 
ones, is given by lSchulten fc Gordon! ( 119751 ). 



The semiclassical approximation for the square of the 6-j symbol is 



i' £ L 
j J J 



-\ 2 



,V12ttV 



cos 



(e + 7r/4) 



1 



2AttV 



;i + sin2e) 



(16) 



and zero when < (IPonzano &: Reggd Il968l ). The semiclassical phase G is defined by 
i^kiikdik where jik is the length of the edge opposed to edge ik and 6^, is the dihedral 
angle between faces intersecting on the edge opposed to ik. Wigner's estimate is obtained 
when the oscillatory factor, from the above formula, is neglected. 



Using notation z = L/n, cos 771 
is calculated as 



n and cos 772 = i'/n, the volume of the tetrahedron 

1/2 



V = { [sin^(77i + V2) ~ z] [z- sin2(?7i - 772)] } 
The generalized character function Hji is a solution of a differential equation, 



;i7) 



dHjL{oj) 
2cotcu + 



4j(4j + l) 



L{L + 1) 
sin^ u 



derived from the corresponding differential equation for Gegenbauer polynomials. A more 
convenient form for Equation (fT8|) is obtained by setting / = sinu HjL{uj) and u = 7r/2 + x. 
The resulting equation 



fix) + 



(2j + 1)' - 



L{L + l) 

COS^ X 



fix) = 



(19) 



is a one-dimensional Schrodinger equation for a particle moving between — 7r/2 < x < 7i/2 
in a symmetric potential. A simple WKB solution for this equation can be constructed to 
obtain the approximation 



hUx) 



2 sinx 



[sin^ -(L/n)^] 



-1/2 



(20) 
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Finally, by using the approximation Equation f|T7|l for the 6-j symbol and Equation fl20|) 
for the generalized character in summation Equation f|T^ , a semiclassical approximation for 
the angular momentum changing probability is obtained as 



vrn^^ sm x 



j dz [[z - sin^(?7i - 772)) (sin^(?7i + 772) - z) (sin^ x - -2)] 



-1/2 



P*^(£/n,£'/n,x) 

TT'/i" sm Y / ^ ■ 

(21) 

T his is precisely the sarn e formu la obtained from purely classical phase space arguments 



m 



Vrinceanu fc Flanneryl (j2001al ). The integration is limited to a proper domain, where 



the argument of the square root function is positive, otherwise, in the classical limit, the 
probability for transition is zero. Moreover, the integral can be calculated in terms of the 
compl ete elliptic integral K{m) = Jq^'^{1 — msin^x)^^/^ dx (formula 3.131.4 of Gradshteyn 
(I2OOOI )). to finally obtain the semiclassical approximation Equation ([6]). 
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Fig. 1. — Probabilities for a dipole allowed transition in the n = 40 degenerate hydrogen 
manifold, for a 36 — )■ 35 transition, as a function of the scaled impact parameter b/an, for a 
fixed projectile velocity f = 0.1 in atomic units, a„ = n'^aQ. Exact quantum probability (red), 
semiclassical (blue), and Born approximations (black) are compared. Th e dotted line marks 
the po sition of the inner cut-off radius used in the Born approximation in lPengelly &: Seaton 
fll964h . 
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Fig. 2. — Integral factor Ii^f^, in Equation (5) calculated rigorously and using a simpli- 
fied semiclassical expression, as a function of the scaled angular momentum transfer. The 
continuous line, the empty circles, and the full circles are obtained from the expressions in 
Equations (7), (6), and (2), respectively. 
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Fig. 3. — CTMC calculation of angular momentum mixing rates scaled by 'on^'KMj {2kBT), 
as a function of the scaled final angular momentum. Dots, squares, and diamonds are the 
results for the three cases discussed in the text, with their statistical errors. Solid lines are 
predictions given by Equation (1131) . integrated over the same final angular momentum bins 
used in simulations. 



